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Theory of strong outcoupling from Bose-Einstein condensates 
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We study the dynamics of a magnetically trapped Bose-Einstein condensate in the presence of 
an external electromagnetic field coupling trapped, untrapped and antitrapped Zeeman sublevels. 
For large condensates an approximate analytical solution of the coupled Gross-Pitaevskii equations 
is given in the regime of strong outcoupling. The theory is developed for the cases of rf-outcoupling 
within a hyperfine manifold of states, micowave-outcoupling connecting states in two different hy- 
0^ ' perfine manifolds, and Raman outcoupling. 
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I. INTRODUCTION 



One of the most promising future applications of the Bose-Einstein condensates of magnetically trapped atoms is 

the possible realization of efficient atom lasers, whose experimental proof of principle has already been given ||l|. An 

^^ ' essential element of any atom laser is the process of outcoupling from the trapped Bose-Einstein condensate. In H] 

an r/-field was used to coherently transfer atoms from the trapped state to untrapped magnetic sublevels within the 

-vj ' same hyperfine manifold of states, enabling them to leave the trap under the influence of gravity. Another outcoupling 

^ , scheme presently being explored experimentally H is based on Raman transitions in an external optical field. 

(M^ The theory of outcoupling from Bose-Einstein condensates has been developed in a number of papers [p|,HJ§-@l ■ The 

I/"} _ theory is based on the Gross-Pitaevskii equations for the macroscopic wave functions of the atoms in the trapped 

and untrapped magnetic sublevels, coherently coupled by the externally imposed electromagnetic field. In |^ the 

coupled Gross-Pitaevskii equations for two magnetic sublevels were solved numerically in 1 dimension for a number 

of cases, and various dynamical regimes were identified. In [J] a detailed comparison of numerical and experimental 

results was given. In [pj the regime of weak outcoupling was considered and in this limit an approximate analytical 

solution of essentially the same theoretical model as in g] and Q was given. The weak outcoupling regime is defined 

C^ ' physically by the condition that the Rabi oscillations induced by the electromagnetic fields in a certain resonance 

Ch \ zone within the condensate are slow on the time scale on which the driven atoms leave that zone. An experimental 

I ■ study P] of this regime was recently performed. In this limit the outcoupling proceeds at an approximately time- 

'^ independent perturbatively calculable rate ||] , and the reduced effective condensate dynamics can be described by a 

C , Markoff process, similar to the familiar Markovian description of a lossy laser mode. The opposite regime of strong 

outcoupling, where the atoms move very slowly on the time-scale of the Rabi cycles, has also been seen in numerical 

simulations y,H but has not yet been studied theoretically in detail. In fact, the experimental realization of an atom 

laser in m operated in this regime, which is therefore of practical as well as theoretical interest. In the present paper we 

give an approximate analytical solution of the coupled Gross-Pitaevskii equation in the regime of strong outcoupling. 

This will be done by a method similar to that of iiq], |1^ for the solution of the time-dependent Gross-Pitaevskii 



C^ ' equation. 



II. MODEL OF THE OUTCOUPLING PROCESS 



The model of the outcoupling process we shall analyse in this paper is essentially that of ref . |^ , but with allowance 
for the action of gravity, which is actually very important [M , at least in the case of microwave or r/-field outcoupling, 
for extracting the outcoupled atoms from the trap. Thus we consider a coupled set of Gross-Pitaevskii equations 
for the macroscopic wave function -(/;„ in the different internal atomic states, e.g. Zeeman sublevels, labelled by m. 
We shall assume that the interaction energy between the atoms is a functional of the total density of the atoms 
IV'P = X)m iV'mP only and not dependent on the internal state m. This assumption is rather well satisfied within 
a hyperfine manifold of fixed total spin quantum number F , but we shall use it as a simple model assumption also 
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for levels with different F-values. Then the coupled Gross-Pitaevskii equations for r/-outcoupling from an F = 1 
hyperfine manifold or microwave oucoupling between two different hyperfine states take the form: 

ini;,n = (-^^ + Vmi^) + UoHA ^m + ?ir!*e*'^* ^m+i + nne-'^'^/j„,-i (2.1) 

M is the atomic mass, Uq — '^^^j " describes the effective low-energy scattering potential with scattering length a. We 
shall restrict our discussion to the case of repulsive interaction a > 0. 

For r/-outcoupling from an F = 1, mp — —1 condensate, like the ^^Na condensate in ||l| (cf. fig|^(a)), m — —1, 0, 1 
labels the 3 magnetic sublevels of the F = 1 state with trapping potentials 

V±i = ± (v{0) + Y'^^rA + Mg(z - zq) (2.2) 

Vo = Mg{z - zo) . 

Here and in the following we neglect quadratic Zeeman shifts. Thus m = and m — 1 refer to the untrappcd and 
antitrapped magnetic sublevels, respectively. Here we included the gravitational potential and denote by V{0) the 
offset of the trapping potential for the to = — 1 condensate in the geometrical center of the trap at r = 0. For 
simplicity we shall assume that the magnetic trapping potential is isotropic. The r/-frequency is denoted by uj and 
ri = guBlBrfl/V^fi is the Rabi frequency at resonance due to the magnetic r/-field, which is spatially independent 
over the condensate. 

Similarly, for r/-outcoupling from an F = 2, mp = +2 condensate (cf. fig.nVb)), like the ^"^ Rb condensate, 
771 — —2, —1, 0, 1, 2 labels the 5 magnetic sublevels of the F ~ 2 state with trapping potentials 

V±2 = ± (v{0) + Y^^rA + Mg(z - zq) 

V±i = ±i (v{0) + Y^^rA + Mg{z - zo) (2.3) 

Vo = Mg{z ~ zo) 
The coupling term between magnetic sublevels in eq.(2T) now reads |1J], |13|| 



i^Q(v/(2 - TO)(2 + TO + l)e*'"Vm+i + V(2 + m)(2 - to + l)e-^"Vm-i) 

with H = g[iB\Brf\l'h. 

For microwave outcoupling between two different hyp erfine states it is sufficient, due to energy conservation, to 
consider only two components tti = 1, tti = 2 in eq. ( p.l| ) which refer to the trapped and untrapped (or antitrapped) 
states, respectively. With the magnetic dipole matrix-element /ii2 between the two states , the Rabi frequency is 
^ = Mi2|^|/2?i where the magnetic microwave field is again spatially independent across the condensate. In this case 
the potentials Vm are 

Fi = yGJ't' + 1/(0) + Mg{z - Zo) 

where p = and p > for outcoupling to an untrapped state and an antitrapped state, respectively. 

For Raman outcoupling (cf. fig.y(c)) ttt. = 1, to = 2 again refer to the trapped and untrapped (or antitrapped) state 
respectively, and fl is replaced by H, ■ e^^^ ^ the space-dependent effective Rabi frequency of the Raman transition, with 
net energy and momentum transfer to the atom of Tilj and ?ifc, respectively. In a three-level approximation for the 
Raman transition fJ — —^j-^- Here ili, 0.2 are the Rabi frequencies of the two lasers with frequencies wi, L02 driving 
transitions between levels 1 and 3, and levels 2 and 3, respectively, with to = toi — 0J2 and k — ki — k2. fl depends 
on the detuning S = {E3 — Ei)/fi — uji between the level-spacing F3 — Ei and hui, which is a ssumed to be large to 



avoid the population of the auxiliary level 3. The potentials Vm in this case are as in eqs. (2.4). 

In the following we shall assume that the electromagnetic outcoupling fields are switched on at time t — 0. We shall 
consider the case of sudden switch-on in a time interval short compared to the final Rabi period. Before the switch- 
on the condensate state is assumed to be in equilibrium and described to sufficient accuracy by the Thomas-Fermi 
approximation, while the outcoupled states are unpopulated. 



III. RABI-OSCILLATIONS IN THE CONDENSATE 

For sufficiently strong external electromagnetic outcoupling fields the Rabi oscillations between the coupled internal 
atomic states become so fast, that the center of mass motion of the atoms can no longer follow. In this r egim e a 
combined Thomas-Fermi and Raman-Nath approximation applies, in which the kinetic energy terms in eq. ( ^.l[ ) are 
negligible. The Rabi oscillations within a given F hyperfine manifold can then be solved analytically. We consider 
the cases F — 1 and i^ = 2 in turn: 

1. Rabi oscillations in an F = 1 manifold: 
Before starting let us shift the origin of the 2-axis into the minimum of the total external potential (trap and gravity) 
seen by the condensate in the ?7i = — 1 state via the replacement 



z^ S^ z + g/u,^ . 




Furthermore let us choose the zero of the gravitational potential in z = by adopting zq = 


—g/uj^. With the notation 


f 2 = a;2 _^ y2 _^ ~2 ^ -(>(Q) ^ y,^^-; ^ 0) = 1/(0) + Mg^/2uj^ 


(3.1) 


we then have 




V^, = ^cD'r' + ViO) 




K, = Mgz 


(3.2) 


Vi = u!^f^ - V{Q) + 2MgS 





(seefigl). 

Furthermore let us proceed to an interaction picture by splitting off the main frequencies according to 

^_i = e-*(A'+V'(o))tA^_^ 

iPo = e-*(^+^(''»*/^e*'^*V;o (3.3) 

Initial conditions have to be specified next. We shall assume the electromagnetic outcoupling fields are switched on 
at t = in a time interval short compared to the final Rabi period. The initial condition at t = in the Thomas-Fermi 
approximation Il4|, |l5[ then is 

V^o = = Vi (3.4) 

\ Uq 2 

where the chemical potential /i is determined by the number of atoms in the condensate via 

N^ ^i2^^f/y [Mcor^^ (3.5) 

15C/o 

Neglecting the kinetic energy terms in a zeroth-order approximation, the coupledGross-Pitaevskii equations reduce 
to a Rabi level-flopping problem, in which the local quantity \ip{f, i)p is constant in time. In the next (first-order) ap- 
proximation \tjj{f,t)\'^ will be allowed to change slowly in space and time and will then be determined self-consistently. 
Therefore, in the present zeroth order approximation we have to treat |?/'('?,t)P as constant in time, but arbitrary. In 
a vector and matrix notation, with 



the coupled Gross-Pitaevskii equations to zeroth order then reduce to the simple form 




e n 

iip^ { n e + A n 

n e + 2A 



^ 



(3.6) 



with 



neif,t)^Uomf,t)\'-f,+ fu>'P^ 



hA{r) 



M 



-2-2 

u! r 



MgS - V{0) + huj 



(3.7) 



Its solution is recorded in the appendix, together with the corresponding solutions for the other cases considered 
below, in particular the case F = 2. Let us define the regime of strong outcoupling by the condition 



ri>|A(r) 



(3.8) 



throughout the condensate. If the r/-frequency lo is chosen in such a way that the resonance condition A(r) — 
is satisfied somewhere inside the condensate, a sufficient condition for (3.8) to hold is f2 ::^ /i. Since in the regime 
of strong outcoupling the detuning A is negligible compared to 17, the solution of eq. (3.6) becomes very simple and 
reads 



^(f,t) 



VW^e 



^iJ*dr[e{f,T)+Air)] 



1 + cos V2nt \ 

-~iV2sinV2nt 

-1 + cos Vint I 



(3.9) 



Fig.0 shows the probabilities Itpmirji)]"^ /\ip(r,t)\'^ as a function of time. Thus the Rabi oscillations in the strong 
outcoupling regime occur with approximately constant amplitude throughout the condensate, but with a space- 
dependent phase. In the next order approximation the gradient of this phase gives rise to a particle current, which in 



turn changes the number density |^(f,t)p by particle number conservation. This will be considered in section [V, 
2. Rabi oscillations in an i^ = 2 manifold 

Let us assume the condensate is in the F ^ 2, nip ~ 2 state. We perform an analogous sequence of steps as in the 
prev ious subsection for the case F = 1 Introducing z = z + qjuP' and choosing zo = —g/ui"^ we have in the notation 

(PD 



V2 = ^LO^f^ + f (0) 

1 fM , , ~ , 
Vq = Mg~z 



V-i 



Lu^r^ + V{0) 



-MgS 



-Mgz 



(3.10) 



M 



V-2^-[ Y^^r^ + ^(0) ) + 2M55 , 



Let us split off the main frequencies according to 



^^ ^ ^-i{^L+ViO))t/h^~^{m^2)uJt^^ 



(3.11) 



Using the fact that |^(i)P = J2m. \4'm{t)\'^ is locally conserved in the combined Thomas-Fermi and Raman-Nath 
approximation and left arbitrary at this stage (to be determined self-consistently later on) we obtain the coupled set 
of equations, in matrix notation with 



i' 



W-2/ 



iip 



Is 


n 








\ 


n 


£ + A 


v/f" 











V 2 


e + 2A 


V 2 











Y 2 


e + 3A 


n 


\o 








n 


e + 4A/ 



V' 



(3.12) 



where 



fe = t/o|Vp-/i+f<I'V~' 



aA 



-Uf^V + no) 



—Mqz + ?icj . 
2 



(3.13) 



The solution of ( 3.12 ) is given in the appendix. In the hmit of strong outcoupHng it takes the form 



V' = 2^^^^ 



-if dT[s(f,r)+2A{f)] 



( cos f2i + i cos 2^t +1 \ 
—i sin nt — ^ sin 2f7i 
^(cos2m-l) 
isin J7i — ^ sin2f7i 

\ - cosnt + \ cos 2r2i + 1 y 



(3.14) 



Fig.H shows the probabihties |V'm(f, t)p/|-0(r, t)p as a function of time for this case. 
3. Rabi oscillations for microwave and Raman outcoupling 

For these outcoupling mechanisms the internal state of the condensate atoms is coupled to a single magnetic sublevel 
in another hyperfine manifold and a simple two-level Rabi oscillation results. As before we introduce z = z + g/uP' 
and pick zq = —gjuP'. Then 



v^i = |:^v+no) 



M 



V2^-p[ Y^^r^ + V{0) ) + Mg{l + p)z . 



(3.15) 



We define V'm via 



^1 



-i{p+v{o))t/h 



^1 



lP2 = g'i(M+V(0))t/h+iujt 



i'2e 



ikz 



(3.16) 



where k ^ Q only in the case of Raman outcoupling. Then the equations in combined Thomas-Fermi and Raman-Nath 



approximation in a matrix representation with ip = 



'W2^ 



read 



iip 



e n 
fl e + A 



^ 



(3.17) 



with 



M 



2 1,2 



UA 



h'k 

2M 



(1 + P) ( Y^"^^"^ ~ ^^^^ ^ ^^^"^ ]+'hLo. 



(3.18) 



Their solution in the limit of strong outcoupling fJ 3> A is 



^(i) = ^ww 



cos57t 
—i sinf7f 



-i/Jdr(e(r,r) + iA(r)) 



(3.19) 



IV. CURRENT OF ATOMS WITH OUTCOUPLING FIELDS ON 



In the preceding section the kinetic energy term in the Gross-Pitaevskii equation ( ^.l| ) was neglected. No spatial 
transport of particles can occur in this approximation, the atoms are assumed to be at rest. The most important 
consequence of the kinetic energy term in eq. (2.1) is that a spatial transport of the atoms occurs. It obeys the 
conservation law 



dt 



V • j = 



with the current density 



3 = 



h 



2Mi ^^ 



{rra^^m - i^mVrra) 



(4.1) 



(4.2) 



We shall now calculate this current using the results obtained in the preceding section. 

The present discussion applies to all cases 1) to 3) of sec. III. Let us begin with the case of r/-outcoupling from a 



F = 1, mp = — 1 condensate. The state- vector ip is space-dependent due to its proportionality to ^|i/)(f,t)p and its 
dependence on e and A which both depend on f . As shown in the appendix ij} can be written as 

i! = \il}\e'^'p (4.3) 

where the vector Lp is space-dependent only via its dependence on the space-dependent detuning A(r) and has unit 



norm. The phase 4>^ according to eq. (3.9), satisfies in the local rest frame of the atoms 



^ = -(£(f,i) + A(f)) 



Inserting (4.3) in the definition of the current density results in the expression 

j{f,t) = v{f,t)\ij{f,t)\^ 
with the local velocity field in the laboratory frame 



v{f,t) 



M 



2Mi 



Yl ifrn'^fm " 'Pm^'P*m 



As is shown in the appendix the second term on the right hand side of eq. (4.6) is given by 

-sin{V2rtt){uJr ■ 



sin{V2rtt)(u}r — — ( 
2s2 LU 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



For strong outcoupling it is rapidly oscillating on the time-scale tD~^ of the center of mass motion and can therefore 
be neglected in the following. If taken into account it would describe a small oscillatory component in the temporal 
evolution of the center of mass and the radius of the condensate. 

The equation of motion of <f) in the laboratory frame includes the kinetic energy Mv'^/2 = ?i^(V(^)^/2M in addition 
to A and e and reads 



^ = - (h{V^r/2M + e(f , t) + A(f )) 



dt 
It has to be solved together with the continuity equation 

Before doing this it is convenient to make a transformation to the rest-frame of the center of mass 

■il;{r,t) = e " ^ Jo ^ ' ^ h y ' ^^[x,t) 



(4.8) 



(4.9) 



(4.10) 



with ^px — I'ipxle^'^" and 

■ * 

X = f — f (t) . 

The continuity and phase equations in the center of mass frame read 



dt 



V ■u\^:^{x,t)\' = Q 



with 



and 



«= J^^<f>x{x,t) 



(4.11) 



(4.12) 



(4.13) 



dt 2M 



iV4>^f + Mf (t) ■x-eif*{t) + x,t)- A{f*it)+x) 



With the center of mass motion f* — —^gt^Bz and eqs.(3.7) the phase equation simphfies to 



(4.14) 



The ansatz 



dt 2Ar ^ ' h ' ' h h 



IV-.P = {A{t) - B{t)x') e {A{t) - B{t)x') 
0, == a{t) - I3{t)x^ 



solves equations (4.12), (4.13), (4.15) exactly if the coefhcients satisfy the differential equations 



M ^ 

a = A 

n 



and, by the normalization condition N = J cfix\tp: 



2 

X I 5 



2/5 



The initial conditions follow from the Thomas- Fermi wave function at t = 0. They are 

5(0) = ^, «(0)=0 = /3(0). 



The time-dependent Thomas-Fermi radius is related to B{t) via 



TTpit) = rTF(O) 



m 

B{t) 



1/5 



Its equation of motion follows as 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



fTF{t) 



4fW 



= 0. 



It can be integrated once to give 



'TF 



W + ^ 



2(1-2 ^5,^(0) _ 2 
3~ rlp{t) ~ 3 



o'^''^tf(0) . 



(4.21) 



(4.22) 



At large times, where rxpit) ^ »^tf(0), the Thomas-Fermi radius expands balhstically at a constant rate 

r-TFioo) = J^CJrTFiO). (4.23) 

At small times it expands with acceleration uj'^rTFiO) 

1 



rTFit)^rTF{0)il + ^Cuh^j. (4.24) 

The Thomas-Fermi radius obtained by integrating eg. ( 4.22 ) is plotted in fig.|5| as a function of time. 
The result for the condensate density in the laboratory frame is 

mt)\' = -^- ^^^' ,f ^ e 4^(t) - r + -gee, . (4.25) 



The phase of the macroscopic state vector ?/> can also be expressed by r^p (t) as 



m-y{^)-^-^i-'^i~^~^T drr^^yp--X\dr^'^Ar-rii)f^y-J^^ (4.26) 



^Jo V4f(^) / 2^^ "•" dt TTFiO) 

Thus under the influence of the strong electromagnetic outcoupling fields, switched on at i = 0, the whole condensate 
undergoes rapid Rabi oscillations through all the coupled states while freely dropping out of the trap by gravity, and 
at the same time expanding due to the repulsive interaction. Switching on the external electromagnetic field therefore 
effectively switches off the trapping potential, up to the oscillatory term (4.7) whose influence is small if cj << ft, 
because it averages out in the Rabi oscillations through trapped, untrapped and antitrapped states of a hyperfine 
manifold. 

The results given so far have been written down for r/-outcoupling from an _F = 1, niF = — 1 condensate (OI 1), 
but are easily transferred to the other cases considered in section (UT): 

RF-Outcoupling from the F — 2, rap = 2 state: 

The velocity field v is still given by eq.(4.6), but the second term on the right hand side of that equation is changed 
and becomes now 

- ^sm(rit)((I'f - l-e^). (4.27) 

It is again rapidly oscillating on the relevant time-scale uj and hence negligible. The equation of motion of the phase 
becomes 

After inserting e and A from eq. ( 3.13| ) , the phase equation reduces to the same form as in the case of the F = 1 
manifold and therefore the results for f*{t),rTF{t), IfAP and 4> are unchanged. 

Microwave and Raman outcoupling: 
We give the results for Raman outcoupling. To obtain the results for microwave outc oupl ing from the following 
formulae on has to set A: = there. The velocity is again given by the general expression (4^), with the second term 
now of the form 

(1 -t- p)uj , „ , , Q 

^— ^sz7i 2m cDr - ^e, , 
41 Z uj 

which is again rapidly oscillating on the time-scale Hi and therefore negligible. The phase equation then takes the 
form 

with 



he 



nuj/2 + nuo\i'\^ - ti- 



h -(1 -p) — OJ r^ 

AM 2 ^ ^'2 



+ i(l+p)M.gz-i(l+p)f(0), 



(4.30) 



Let us first consider the case p = 1. There the excited state is antitrapped and has a magnetic moment precisely 
opposite to the trapped condensate state. The trap potential thus averages out in the Rabi oscillations. We transform 
to the center of mass frame with 



hk 

M 



gte-z 



(4.31) 



There we can again make the ansatz (4.16) and solve for A{t), B{t) as before. Finally, transforming back to laboratory 
coordinates we obtain 



mt)? 



15N r^pit) - (r + ^e.g^ - ^t) 



Stt 



' TF 



it) 



(4.32) 



with the previous result ( 4.21 , 4.22 ) for the equation of motion of rxp (t) ■ 

In the case p ^ 1, to which we now turn, the trap potential does not average out in the Rabi oscillations and 
contributes to the center of mass motion and the spreading of the condensate. Now the transformation (written for 
p < 1) f = f*{t) + X with 



1 + P gez 
1 — p (D^ 



r-(t) = -i^^ 1 



1-p 



Lot 



hk 

M 



-uit 



(4.33) 



introduces the center of mass rest-frame. The phase equation in this frame takes the form 



dt 



2M 



i^k? 



Uo 



H.\ 



l-p M_2 2 

— — —LU X 

2h 2 



l"pM_2-,2 1+P 



fi hk^ 1 + p - 



2h 2 



2h ^ h AM 2n ^ ' 2 



(4.34) 



Again the continuity equation and the phase equation can be solved by an ansatz of the form (4.16). The solution for 
the condensate number density with the required initial conditions takes the form, in the laboratory frame. 



W{tt 



IbN rlp{t)-{f-f*{t)y 



Stt 



4fW 



e((f-P(t))^-4^(i)) 



where the Thomas-Fermi radius rTP [t) satisfies 



fTF{t) - 0? 



4f(0) \-p 



'TF 



(t) 



TTpit) 



(4.35) 



(4.36) 



It is not influenced by the center of mass motion f*{t). The phase of the state vector -0 expressed in terms of rxpit) 
is 

hk'^ l-fp ~ , , uj\ 
-4M+^^(0)--h 



M /■* 



drii^-l]+^ir-rit)f^ln 



TF 



(r) 



M 
2h 



dt'" rTF{0) 



M ^ , , - 
+—r*{t) r 
n 



(4.37) 



where energy conservation of the center of mass motion led to the cancellation of some terms. Thus, for p < 1 the 
condensate remains bound, even in the presence of the electromagnetic outcoupling field. Its center of mass oscillates 



with the frequency \ —^'j^ around the new equihbriuni position at z = — j^ ^- ^^'^ Thomas-Fermi radius also 



oscillates, for small amplitudes with frequency \/f (1 ^ p)^ around its new equilibrium value (see figra) 



2 ^^/^ 



rTF=rTF(0)( Y^j ■ (4.38) 

For p > 1 the condensate becomes antitrapped by the outcoupling field and is driven apart by the magnetic field of 



the trap. Its center of mass trajectory can be inferred from (4.33) by analytic continuation to p < 1 as 



sinh ( \/^^-T^ujt 



Its Thomas-Fermi radius (see fig.0) grows with acceleration ^^uj'^rTFiO) for small times cut <$:. 1 

TTFit) = rTF(O) n + ^^^LdhA , (4.40) 

and exponentially at large times 



TTFit) ~ exp ( \l^^u;t J (4.41) 

provided the condensate has not yet spread beyond the parabolic part of the trapping potential. 

V. OUTPUT 

We assume the outcoupling field is turned off rapidly when, in a Rabi cycle nearly all atoms are in the antitrapped 
state. To determine their macroscopic wave function afterwards we have to solve their Gross-Pitaevskii equation 
again. 

1. F = 1 manifold 
For outcoupling from the _F = 1 manifold 

iHi = ( ~^^' ~ y^'^' + 2M5f - /i - 2y(0) + 2hu; + C/olV^iH i>i (5.1) 

where atoms remaining in the to = —1,0 states have been neglected in the total density. This is justified if the 
outcoupling field is switched off at a time T satisfying 

V2nT ^ {2n + 1)tt (5.2) 



with integer n. To solve ( |5.l[ ) we can use the same method as in section IV. The center of mass r*(i) for t > T 

satisfies 

.. ^ 

r = wf* — 2gez . (5.3) 

With the apropriate initial conditions at time i = T it is given by 

f*{t) = ^^ (l " (l + ^^-^) cosh(I'(t - T) - — sinh(I>(i - T)) , (5.4) 

i.e. the combined action of gravity and the antitrapping potential pushes the center of mass of the outcoupled 
condensate vertically downwards. Meanwhile its Thomas- Fermi radius keeps expanding according to the equation of 
motion 

f^HO--^(^+rTF(t))=0 (5.5) 
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which corresponds to eq. (4.36) with p = 3. It has to be solved with initial conditions at time t — T, ttf — ttf{T) 
and 



1/2 

rTF{T) = \l ^OjrTF{0) ( 1 - ^^^9;^ ) (5.6) 



i-M^-m 



following from the solution for t < T of eq. ( 4.21 ). The solution for rTF{T) is plotted in fig.|§. For times ujt ^ 1 



but still small enough so that the atomic cloud has not left the parabolic part of the trap, the cloud's center of mass 
coordinate and its radius both grow e xpon enti ally with the sam e rat e uj. The density and phase of the outcoupled 



wave function are then given by eqs. (4.35), (4.37), where eq. (4.37) has to be taken for p = 3. Fig.y gives a plot 



of the number density \ip\ in the antitrapped stated. Plotted is |V^ in units of 157V/87rr|,p.(0) as a function of the 
radial coordinate in the center of mass frame in units of rTF(O) and f2i. In order to portray the influence of the 
Rabi-oscillations and the motion in the trap potential together the ratio Jl/w in this and in the following two plots 
has not been taken as large as would be desirable for the clear separation of time-scales we have assumed. In this plot 
cS/ri = 0.1, and the duration of the radio-frequency pulse is uiT = 2. 

2. F = 2 manifold 
In the limit A/i7 -^ a complete transfer of the occupation probability from the trapped mp = 2 state to a single 
antitrapped state (cf. fig.^) occurs at the times flT = {2n + 1)tt when all atoms are in the antitrapped niF — —2 
state. We can take over the preceding analysis with minor modifications. The Gross-Pitaevskii equation of the atoms 
in the uif = — 2 state 



inij_^ = ( -2^V2 - —6^2-2 ^ 2MgS - /i - 2T/(0) + AHu; + l/o|V^-2p ) 4'-2 (5.7) 



takes the same form as (|5.lD if there the replacement huj — > 2huj is made. The result for the density in the tuf = —2 
state is plotted in fig.lO. 

3. Microwave and Raman outcoupling 

In the limit A ^ J7 all atoms in the condensate are in the antitrapped state m = 2 (for p > 0) at times ilT = '^"^^ tt. 
The subsequent evolution of ip2 is governed by 



(5.8) 
-M+^ + ?^^-(l+p)n0)jvi2. 

The center of mass satisfies the equation of motion 

r = pui'^f* — (1 + p)gez (5.9) 



with initial conditions following from cq.(4.33). The Thomas-Fermi radius satisfies 

TTFit) - CD' (^^^ + pTTFit)^ = (5.10) 

with initial conditions rxpiT) and 



3 I r^FiT)) 2 ^v4f(0) 



rrFiT) = urMOhro 1 - ^F^ + ^ T^ 1 (5.11) 



following from eq. ([4.36|) . The solution for p = 1 is shown in fig.|8|. The density |V^2p is then given in terms of f*{t) 
and rTF(t) by eq. ( [4. 351 ). For the purpose of generating an output pulse with a small Thomas- Fermi radius the case 
p = is of particular advantage. The number density in the excited state for this case is plotted in fig.y_l[ In thi s plo t 



it is assumed that il/uj = 5.91 and flT = 157r/2. Finally, the phase is obtained as an expression similar to eq. (4.37) 
namely 
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{f,t) 



+ {i+p)v{Q)-nLo^ ' 



^ ./o 



2Af 

t 



m 



'TF 



(r) 



dr (l|:n^-l)+i^(r-r*(i))2^1n 



M 
2h 






(5.12) 



r*(t) ■ r 

h 



where rj.p{T) in the time integral has to be taken as solution of eq. (pjq) for < r < T and as solution of eq. (5.1C) 
ioiT <T <t. 



VI. CONCLUSION 



In this work we have studied outcoupling from magnetically trapped Bose-Einstein condensates via electromagnetic 
fields resonantly coupling either all the magnetic sublevels in the hyperfine manifold of the trapped state or the trapped 
state and a state in a different hyperfine manifold. Simplifying assumptions we made were (i) isotropic trap potentials, 
(ii) equally spaced magnetic sublevels, i.e. neglect of quadratic Zeeman shifts, (iii) interaction of the particles via 
their total density, (iv) the validity of a dynamical version of the Thomas-Fermi approximation for the condensate, 
in which the kinetic energy per atom is small compared to the chemical potential and (for h — 1) the Rabi- frequency, 
and (v) the negligibility of thermal excitations. With these assumptions we have given an analytical treatment of the 
problem of strong outcoupling, where everywhere in the condensate the Rabi frequency dominates over the space- 
dependent detuning, which is at most of the order of the chemical potential. The equations of motion derived in 
that hmit (eqs.( [4.12 ), ( [1.13D , and ( |4.15| ) for rf-outcouphng, eqs.( [4.12D , ( f4.13| ), ( f4.29| ), ( |4.30| ) for microwave and Raman 
oucoupling) could be solved exactly. For its quantitative validity our analytical treatment depends crucially on the 
assumptions (ii), (iii), (iv), and (v), but we expect that qualitatively our results should even be applicable if these 
assumptions are not well satisfied. On the other hand the method we use can be generalized to anisotropic harmonic 
potentials, so that assumption (i) could be relaxed, but at the cost of more tedious analytical expressions. Our method 
of solution expresses the densities and the phases of the time-dependent macroscopic wavefunctions in the coupled 
states in terms of the center of mass coordinate and the Thomas-Fermi radius of the condensate and theoutcoupled 
state and provides the ordinary second order differential equations satisfied by these quantities. 

The center of mass motion is in all cases very simple. It is, of course, independent of the interaction between 
the atoms and driven by the momentum transfer from the outcoupling field to the atoms in the case of the Raman 
outcoupling, by gravity, and by the time-averaged trapping potential seen by the atoms. The latter vanishes for 
resonant rf-outcoupling within the hyperfine manifold of the condensate, i.e. the center of mass in this case is freely 
falling after the outcoupling field is switched on. If the outcoupling field is off-resonance, and in the cases of microwave 
or Raman outcoupling, the time-averaged trapping potential is, in general, non-zero (cf. eq.( 4.33| ) for the center of 
mass motion in this case) and leads either to a residual harmonic binding of the center of mass to the trap or to 
exponential repulsion from it. 

The second parameter characteri zing t he condensate in the outcoupling fields is the tim e-dep endent Thomas-Fermi 
radius. Its equation of motion (eq.( 4.21 ) while the outcoupling fields are on, eq.(5^) or ( 5.10 ) after the outcoupling 
fields are switched off) are still very simple and can be solved analytically up to quadratures. In figs.H-0 the three 
different types of solutions are portrayed. 

Fig.p shows the expansion of the Thomas- Fer mi rad ius for the resonant rf-outcoupling within the hyperfine manifold 
of the condensate, which is the solution of eq.( 1.21 ). An initial phase of accelerated expansion due t o the mutual 
repulsion of the atoms is followed by a purely ballistic expansion with constant velocity given by eq.(4.23). If the 
outcoupling field is switched off the density of atoms in the antitrapped states enter a phase of expontial expansion. The 
total evolution of rTpit) for this case is shown in fig.H while the complete evolution of the density in the antitrapped 
state is seen in fig.^ for the F = 1 manifold and in fig.|l^ for the F = 2 manifold. The only differencebetween these 
two cases turns out to be the pattern of the Rabi-oscillations, seen directly in figsj|, |j. Fig.0 shows the expansion of 
the Thomas-Fermi radius for a case where the time-averaged potential seen by the atoms is repulsive. The expansion 
again starts out with constant acceleration and crosses over to an expontial phase even while the outcoupling fields 
are on. The third possibility for the evolution of rTF{t) is portrayed in fig.O, which applies to the case where the 
time-averaged trapping potential seen by the atoms in the outcoupling field remains attractive. The Thomas-Fermi 
radius oscillates in this case around its new enlarged equilibrium value either until dissipation takes over (neglected 
throughout this work, which is justified if only a few oscillation cycles are considered) or until the outcoupling fields 
are switched off. If this is done at a moment where most atoms are in the outcoupled state, which may be untrapped 
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or antitrapped, a ballistic or exponential expansion follows the oscillations, respectively. The complete evolution of 
the density of atoms outcoupled to an untrapped state is shown in fig.|ll|. 

It becomes clear from our treatment that external electromagnetic fields inducing rapid Rabi oscillations are excellent 
tools for suddenly changing or effectively turning off the trapping potential, thereby exciting the center of mass and 
the dilation-compression modes of the condensate. A further freedom in the possible changes of the effective trapping 
potential is the introduction of a constant overall detuning of the Rabi oscillations, if lo is chosen in such a way 
that nowhere in the condensate the resonance condition is satisfied. Then only a small fraction of the atoms in the 
condensate is coupled out by each pulse of the external fields. The strong-coupling limit in this case is reached if 
the effective Rabi frequency, including the detuning, is large compared to the variation of the detuning across the 
condensate. Our treatment, in which resonance was assumed, may easily be extended to this case. 

As long as the external electromagnetic field is on, the condensate is not yet coupled out, unless the antitrapped 
states to which it is coupled happens to be more strongly repelled by the magnetic trap than the trapped states are 
attracted, which is the case p > 1 in eq.(|3.15). In all other cases, how many atoms are finally coupled out depends 



crucially on the phase ^T of the Rabi cycle at the moment when the electromagnetic field is switched off. In figs.p[]ll| 
we have chosen the switch-off time T in such a way that the outcoupled state has its maximum occupation probability 
precisely equal to or at least close to 1. A comparison of these figures, and of the formulas to which they correspond, 
shows that outcoupling to an untrapped state as in fig.|ll| rather than an antitrapped state as in figs.y, |lO| is highly 
preferrable to achieve a spatially confined, narrow pulse, as the spatial width of the pulse in the antitrapped states 
is rapidly broadened by the repelling trap potential. However, the coherence of the outcoupled atom wave should 
permit a refocussing of the outcoupled pulse by the techniques developed in atom optics. 

We have followed the dynamics of the atoms only in the region of space where they feel the parabolic part of the 
trap potential. In order to follow them further as they leave the trap and finally become freely propagating atoms it 
is necessary to make detailed assumptions about the trapping potential far from the center of the trap. If this can be 
done for a given experiment one may use semiclassical methods to solve the Schrodinger equation for the atoms [M, 
eventually neglecting their interaction if the atomic beam hasbecome sufficiently dilute. This neglection is possible 
as soon as the rj^^-term in the equation of motion for ttf has become negligible. Our present results for the region 
inside the trapping potential provides the initial or boundary condition which is required for such a further analysis. 
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APPENDIX: 

Here we record the needed analytic expressions for the Rabi oscillations in the various cases considered in this work. 
F=l manifold: 
The ansatz 

^ = e"*/o -^^ [^(*^'")+'^(')l|V;(f,t)|(^(m) (Al) 



with |</?| = 1 transforms eq.(3.6) into 

i(p' = I 1 1 I V? (A2) 




with b = A/il, where ip' is the derivative of f with respect to t = fit. Its solution with initial condition ipm = (5m, i 
is given by 



1 



2 + ^2 



1 + (1 + ^2)cos(V2 + 5^t) + iSV2TS^sin{V2TpT) 

5(1 - cos{Vf+5^T)) - iV2T¥sin{V^TpT) 1 . (A3) 



-1 + cos{V2 + ¥t) 
In the present work we need ip{Qt) for 5 = and 



13 



in order to evaluate the current velocity {h/ M)Im[tp*'V tp]s=o- Using (]A3[) we obtain 



Q 1 / iV2sin{V2nt) \ 

(t|)5=o = 2 M " cos{V2nt) 



35' 



and 



¥—2 manifold: 

The ansatz 



/2 



^^g-/>-[^(^.-)+2A(r)l|^(^^i)|^(f^i) 



with |(/3p = 1 transforms eg. ( 3. 12 ) into 



( -25 10 \ 
1 -(5 a/I 



up 



Vl Vl 

a/I (5 1 



V 



V 1 25/ 



with the same notation as in (A2). The solution with initial condition iy9„i(0) = 5,^,2 is given by 

cos(Ait) — cos(A2t) 



(^(t) = (^(0) + Ci[cos(Ait) + cos(A2r) - 2] + C2- 

+ iC3[Aism(Air) + A2sm(A2T)] + iCa 
with the characteristic dimensionless frequencies 



Aj - A2 



Aism(AiT) — A2sm(A2T) 
Aj - A2 



Ai - 2v/<52 + 1, A2 = v/^2 + 1 



and the coefficient- vectors 



Ci = 



16((52 + 1)2 



C3 



16(52 + 1)2 



/ %5^ + 16(^2 + 5 \ 

' -6(5 ' 

-V6(2(52 _ 1) 

6(5 
-3 



/ 2(5(2(52 + 5) \ 

' 4(52-5 > 

-3^/6(5 

3 

V / 



;C72 = 



1 



16((52 + 1) 



; Ca 



1 



16((52 + 1) 



/ 24^4 _ 9 ^ 
' -65(8^2^3) \ 

3V6(6(52 + 1) 

-30(5 

V 15 ) 



I 6(5(2^2 _ gw 

/ _9(4j2_i) \ 

15V6(5 
-15 

V / 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



The special case 5 = gives Ai = 2, A2 = 1 and leads to eq. (|3.14| ). In the present work we also need (V(/3)5=o and 
hence (-gD^^o- Using (A9)-( A11 ) we obtain 
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and 



Two-level case: 

The ansatz 



leads to 



/ Usin{2nt) + sin{nt)) \ 






iVe 



(cos(20i) - 1) 



-\cos{2nt) + co3{nt) - I 




Im[ip*'Vip]s=o = —sin{Qt)'VA 



^ = e-*/o -^^ ['(^^"^+^(^)/'l|V^(f,i)|<^(m) 



(A12) 



, . _ 1 f \fP^ncos{^\ + 5^ /A t) + i5si n{^\ + 5^ /At) \ 
^^^^ " V¥TA \ -2isin{^\ + 52/4r) ) 



(A13) 



(A14) 



(A15) 



from which we can extract the result for 



/to[</7*V(/7]5=o = — sm(2f)t)VA 



(A16) 
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FIG. 1. Outcoupling schemes for radio-frequency fields (a) micro-wave fields (b) and optical fields (c) 
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FIG. 2. Potentials of the Zeeman sublevels of the F = 1 manifold 
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FIG. 3. Occupation probabilities of the Zeeman sublevels of the F = 1 manifold 



18 



IV, 



ml 












^t 



FIG. 4. Occupation probabilities of the Zeeman sublevels of the F = 2 manifold 
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FIG. 5. The Thomas- Fermi radius rpTit) in units of rpriQ) as a function of ujt during Rabi oscillations within a hyper- 
fine- manifold 
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FIG. 6. The Thomas- Fermi radius rpTit) in units of rpTiQ) as a function of tot during Rabi oscillations between a trapped 
and a weakly anti-trapped state (p=l/2) 
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FIG. 7. The Thomas- Fermi radius rpTit) in units of rpTiQ) as a function of ut during Rabi oscillations between a trapped 
and a strongly anti-trapped state (p=3) 
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FIG. 8. The Thomas- Fermi radius rprit) in units of rFT(O) as a function of uit during {uit < 2) and after (tot > 2) Rabi 
oscillations within a hyperfine manifold or between trapped and antitrapped states with opposite potentials 

FIG. 9. IV'i f for outcoupling from the F — 1 manifold in units of 15N/8Tvr^p{0) as function of the radial coordinate in the 
center of mass frame in units of rTF(O) and time in units Q~^ for the choice of parameters to/Q = 0.1 and flT = 20 

FIG. 10. 1V'-2|^ for outcoupling from the F = 2 manifold plotted as in fig.|for cu/Q. = 0.1 and ilT = 22 

FIG. 11. |^2p for outcoupling from a trapped state to an untrapped state (p = 0) in a different hyperfine manifold plotted 
as in fig.| for ujQ, = 0.169 and J7T = 157r/2 
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